Fano-Kondo effect in side-coupled double quantum dots at finite temperatures and the importance 

of the two-stage Kondo screening 
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We study the zero-bias conductance through the system of two quantum dots, one of which is embedded 
directly between the source and drain electrodes, while the second dot is side-coupled to the first one through a 
tunneling junction. Modeling the system using the two-impurity Anderson model, we compute the temperature- 
dependence of the conductance in various parameter regimes using the numerical renormalization group. We 
consider the non-interacting case, where we study the extent of the departure from the conventional Fano reso- 
nance line shape at finite temperatures, and the case where the embedded and/or the side-coupled quantum dot 
is interacting, where we study the consequences of the coexistence of the Kondo and Fano effects. If the side- 
coupled dot is very weakly interacting, the occupancy changes by two when the on-site energy crosses the Fermi 
level and a Fano-resonance-like shape is observed. If the interaction on the side-coupled dot is sizeable, the oc- 
cupancy changes only by one and a very different line-shape results, which is strongly and characteristically 
temperature dependent. These results suggest an intriguing alternative interpretation of the recent experimental 
results study of the transport properties of the side-coupled double quantum dot [Sasaki et al., Phys. Rev. Lett. 
103, 266806 (2009)]: the observed Fano-like conductance anti-resonance may, in fact, result from the two-stage 
Kondo effect in the regime where the experimental temperature is between the higher and the lower Kondo 
temperature. 

PACS numbers: 72.15.Qm, 75.20.Hr, 73.23.-b, 73.40.Gk, 73.63.Kv 
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I. INTRODUCTION 

In condensed-matter physics, the Fano resonance line 
shape 1 ' 2 is commonly observed in the low-temperature zero- 
bias conductance curves of various mesoscopic and nanoscale 
electronic devices when the energy of a weakly coupled dis- 
crete state is swept across the Fermi level using gate voltages. 
A prototype system where Fano physics may be observed 
consists of a single quantum dot side-coupled to a quantum 
wire^— . The Fano line shape results from an interference of 
the quantum amplitudes for the conduction pathway directly 
through the quantum wire without passing through the quan- 
tum dot and the indirect conduction pathway via the quantum 
dot. The first pathway plays the role of a broad background 
process, while the second corresponds to a resonant scattering 
channel. The conductance as a function of the discrete state 
energy level e is well described by the Fano function 



G(E) 



(E + qY 
'' 1 + E 2 



+ b, 



where 



E = 



r 



(i) 



(2) 



is a dimensionless energy that measures the energy difference 
from a resonance energy e r in units of the resonance half- 
width T, q is the Fano parameter given by the ratio of reso- 
nant and background scattering amplitudes, while a and b are 
some coefficients. 74 This energy dependence holds as long as 
the background conductance is constant over the width of the 
resonance; in experiments, curve fitting with Fano profile is 
typically done in an energy window of the order of several 
times T. 




Figure 1: (Color online) Schematic representation of the double 
quantum dot nanostructure consisting of a quantum dot (d), em- 
bedded between the source and drain electrodes, and a side-coupled 
quantum dot (a). 



The situation becomes more involved in a related structure 
consisting of two quantum dots, the first of which is embed- 
ded in the direct conduction pathway between the source and 
drain electrodes and the second is side-coupled to the first 
one by weak tunnel coupling, as illustrated in Fig. Q] This 
nanostructure, also known as the "side-coupled double quan- 
tum dot", allows to study the interference between the direct 
and indirect conduction pathways (i.e., the Fano effect), as 
well as various correlation effects resulting from the strong 
electron-electron repulsion between the confined electrons, in 
particular the Kondo effect^—. Taking into account only 
the electron levels in close vicinity of the Fermi level, this 
system can be adequately described at low temperatures us- 
ing a two-impurity Anderson model which is a simple ex- 
tension of the single-impurity model that is commonly used 
for studying transport properties of single ultra-small quan- 
tum dots 24 . The side-coupled double quantum dots have been 
intensively studied theoretically 25-33 , uncovering a fascinating 
interplay between many-particle effects and quantum interfer- 
ence, however the first experimental study of this system has 
been performed only very recently (Ref. 341) , Most previous 
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theoretical works consider the zero-bias conductance in the 
zero-temperature limit, while the experiment is decidedly per- 
formed at finite temperatures. It is known that this system 
is characterized by very low energy scales in some parame- 
ter regimes (in particular in the case of the two-stage Kondo 
effect 27 -^5,), thus in order to properly describe its transport 
properties it is imperative to consider thermal effects and to 
calculate the conductance at finite temperatures using a non- 
perturbative method. This is the goal of the present work. 

This paper is structured as follows. In Sec.[II]the model and 
the numerical renormalization group method are presented. 
The thermal effects are first studied in the non-interacting 
model in Sec. [Ill] where we show that at finite temperatures 
the resonance line shape still has the Fano form to a good ap- 
proximation, albeit with temperature-dependent parameters. 
In Sec. [IV] we then focus on the Fano-Kondo effect in the case 
where only the directly embedded quantum dot is interacting 
and experiences the Kondo screening; in this case the reso- 
nance line shape itself is appreciably modified by the com- 
bined Fano interference and Kondo effect. In Sec.[V]we finally 
study the fully interacting case, where completely different 
behavior is found. Finally, in Seel VII the recent experimental 
results are examined and compared to theory, leading to the 
conclusion that the observed line-shape (an anti-resonance) is 
indicative of the occurrence of the two-stage Kondo effect. 



II. MODEL AND METHOD 

The Hamiltonian under study takes the following form: 

H = 5 d (n d - 1) + 5 a (n a - 1) - t ^ (4«a + H.c.) 



+ E e fc c Lc feCT + v J2 (4a + H -' 



(3) 



where n d — d\d a and n a = ^2 a al.a a . The operators d\ 
and a\ are creation operators for an electron with spin a on 
the embedded quantum dot d and on the side-coupled dot a, 
see Fig. Q] The on-site electron-electron repulsion is denoted 
by U a , where a stands for d or a. The on-site energies e a of 
the dots are shifted by U a /2 in order to define the parameters 



$a = £q + U a /2, 



(4) 



which measure the detuning from the particle-hole symmetric 
case which corresponds to S a = S d = 0. The coupling be- 
tween the dots is described by the interdot tunnel coupling t. 
The dot d is assumed to couple to the source and drain elec- 
trodes in a symmetric way, thus it hybridizes with the even- 
parity combination of electrons from both leads 16,36 with the 
strength 



T d = TrpV 



(5) 



which will be assumed to be constant for all energies within 
the conduction band of half-width D (equivalently, the den- 



sity of states in the conduction band, p, is assumed to be con- 
stant, p = 1/2D). Throughout this work, T d will be fixed to 
T d /D = 0.06. 

For the correct description of correlated regimes in clus- 



ters of quantum dot; 



,30.37-51.51-53 



it is necessary to use non- 



perturbative techniques that take properly into account not 
only the local correlation effects on the impurity sites, but also 
the charge fluctuations and inter-site spin-spin coupling in- 
duced by the exchange interaction. One such technique is the 
numerical renormalization group (NRG) 54 69 which involves 
three steps: 1) logarithmic discretization of the conduction 
band, 2) mapping onto a one-dimensional tight-binding chain 
with exponentially decreasing hopping constants, and 3) itera- 
tive diagonalization of the resulting Hamiltonian. The results 
presented in this work have been calculated using the NRG 
discretization parameter A = 4 with N z = 8 equidistant val- 
ues of the twist parameter z using an improved discretization 
scheme^S, with the truncation cutoff set at IOcjat. The ex- 
pectation values are calculated in the standard way, while the 
conductance curves at finite temperatures are computed using 
the Meir-Wingreen formula from the spectral dat a 61,63 ' 71 : 



G(T) = G 7rT d f 

J —i 



did — 



(6) 



where Go = 2e 2 /h is the conductance quantum (i.e., the 
conductance corresponding to full transmission), f(u>) = 
[1 + exp(/3w)] _1 is the Fermi function where (3 = 1/ksT 
and the chemical potential has been fixed at zero energy, and 
A d (uj, T) is the spectral function on the impurity site d. 

It should be noted that there is a significant conceptual dif- 
ference as regards the usage of the Fano formula in the context 
of scattering experiments in atomic physics and the present 
context of nanoscopic transport. In atomic physics, the res- 
onance is observed in the scattering cross-section as a func- 
tion of the energy of the incoming particle while the reso- 
nance energy is fixed. In transport experiments, the resonance 
is observed in the conductance which relates to the scatter- 
ing of electrons at the fixed energy of the Fermi level (this is 
strictly true only at zero temperature; at finite temperatures, 
electrons within an energy window of order fc^T participate 
in the transport) and the independent variable is the energy of 
the side-coupled dot which is swept using gate voltages. In 
spite of the seemingly symmetric roles of parameters e and 
e r in the expressions in Eqs. (HJ and ©, there is a notable 
difference as regards the role of the background channel. The 
background scattering phase shift namely depends on the vari- 
able parameter in the first case, while it is essentially constant 
in the second case (at zero temperature). There may be some 
small variation of the background scattering phase shift when 
the parameters of the side-coupled dot are changed, since the 
two impurities are coupled and there will be hybridization ef- 
fects induced by the side-coupled dot on the embedded dot. 
Nevertheless, this effect is rather small. We therefore con- 
clude that any significant deviation from the universal Fano 
line-shape profile may only result from electron correlation 
effects or from thermal effects. 

Furthermore, it should be remarked that the approach used 
in this work does not take into account that at high tempera- 
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tures the phase coherence necessary for the Fano interference 
may be partially lost, since full quantum coherence is assumed 
in solving the model^^. This implies that in experiments 
further reduction of the resonance amplitude is expected in 
addition to that found in our calculations which only take into 
account the effects due to the thermal broadening of the elec- 
tron distribution function in the source and drain electrodes. 

In the following, we set fee = 1 and the conductance is 
always expressed in units of Go- 



III. NON-INTERACTING CASE 

The Fano formula has been originally formulated for a 
quantum system in its ground state. For a non-interacting 
system (U a = U d = 0), it is thus expected to describe the 
zero-temperature conductance essentially exactly as long as 
the background remains constant, which is the case if the Fano 
resonance is narrower than the Lorentzian spectral peak on the 
embedded quantum dot (of half-width Td). The characteris- 
tic energy scale of the level broadening on the side-coupled 
quantum dot is 



7a(e a ) = 7rA° d (e a )t 2 



(7) 



where A d (u) is the spectral function of the level d in the ab- 
sence of the level a, which is given by 



A» = 



1 



1 



irTd (w - e d ) 2 /T 2 d + 1 



(8) 



in the wide-band limit which applies well here for the chosen 
T d = 0.06D. For later reference, we also give the full spectral 
function of the level d in the presence of the side-coupled dot: 



A d (w) 



1 



n 2 



TrTd [l-(io-e d )(cj-e a )/t^Y + W 



(9) 



with ft = {oj-e a )/{t 2 /T d ). 

The non-interacting case is also a good test for the NRG 
procedure, since a comparison with essentially exact results 
from simple quadrature of the spectral function is possible. 
We find that the results (detailed in the following subsections) 
are correct within the line-width or, more precisely, within 
less than half percent. A similar degree of precision is then 
also expected in the interacting case considered later on, since 
the NRG method functions equally well for any value of U. 



A. Symmetric case, e d = 

Fixing ed = and taking the value of A d (uj) at the Fermi 
level, A° d (0) = l/nT d , we obtain 



r a = 7a (o) = — . 

1 d 



(10) 



The Fano formula applies for r a <C T d , therefore the range 
of validity is expected to be given by t <C T^. In fact, for 



ed = it turns out that the zero-temperature linear conduc- 
tance G(0) = GoTTTdAd(0) reduces exactly to the Fano for- 
mula for any value of the t/Td ratio: 



(ii) 



with E = e a /T a . The asymmetry parameter q is clearly zero, 
thus the Fano resonance takes the form of a symmetric anti- 
resonance. 

The conductance curves for the non-interacting model with 
€d = are shown in Fig. [2^. At all temperatures, the res- 
onance line-shape can be described with excellent accuracy 
using the Fano formula with temperature dependent parame- 
ters which are shown in Fig. |2j). The characteristic tempera- 
ture scale for the temperature-dependence of the conductance 
curves is, as expected, r a . An approximate expression for the 
temperature dependence of the resonance width in the temper- 
ature interval [0 : T a ] is 



r = r a [i + Cl (r/r a ) 2 ] 



21 c 2 



(12) 



with ci = 9.2 and c-2 = 0.46. This temperature-variation rule 
holds only in the t/Td limit: in general, the temperature 
variation also depends on the t/Td ratio. 




Figure 2: (Color online) Non-interacting model with U a = Ud = 0, 
symmetric case, = 0. a) Conductance curves for a range of tem- 
peratures. Fano-resonance-line fits overlap completely with the con- 
ductance curves, b) Fano parameters as a function of the temperature. 



We have verified that for = the Fano curve accurately 
describes the finite-temperature conductance even for t of the 
order of T d (results not shown). 
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B. Asymmetric case, e c i ^ 

For e d ^ 0, the d-level Lorentzian is no longer fixed at the 
Fermi level, which has a number of significant consequences. 
The effective level broadening on the side-coupled dot is re- 
duced and it is given by 



r Q = 7o (o) = 



t z 



i 



r d (e d /r d ) 



l 



(13) 



Furthermore, the background phase shift is no longer tt/2, 
thus the asymmetry parameter q is different from 0. From 
Eq. (O it is clear that the zero-temperature conductance is al- 
ways zero at e a — 0, however the actual resonance energy e r 
is actually shifted away from zero. In fact, within an error of 
order t 2 /T 2 d , the following Fano parameters hold at T = for 
general e d : 

r — r 
q = e d /r d , 



1 



(14) 



a = 



1 + (e d /T d y 
0. 



The results of a temperature-dependent calculation for 
£ d /r d = 0.5 are shown in Fig. [3] We observe a departure 
from the Fano line-shape which grows as the temperature is 
increased; nevertheless, the departure is small and it is still 
meaningful to perform curve fitting using the Fano formula. 
Similar to the symmetric e d = case, the main effect of finite 
temperature is to increase the width of the Fano resonance be- 
yond T a and to reduce the amplitude of the resonance. The 
asymmetry parameter q is, however, nearly constant: it re- 
duces from q = 0.5 at zero temperature to q = 0.47 at 
T = r a , while in the same temperature interval the amplitude 
is reduced by half. The resonance position, as given by e r , is 
also affected rather weakly. Such behavior is characteristic of 
the t <C r d limit. For larger t, the extracted Fano parameters 
q and e r are more strongly temperature dependent, however 
the deviation from the Fano profile remains small (results not 
shown). 

We conclude, fully in line with the expectations, that in the 
non-interacting case the main effect of finite temperature is to 
reduce the amplitude and increase the width of the resonance, 
while the Fano line-shape is largely preserved. As we show 
in the following, this is no longer the case in the presence of 
interactions, except in some special limits. 



IV. INTERACTING CASE WITH U d AND U a = 

The occurrence of the Fano resonance in transport prob- 
lems is always associated with changes of the charge state of 
the weakly coupled discrete state. For this reason, there is a 
significant difference between the U a <C T a and U a ^ r a 
situations. In the first limit, the occupancy changes by 2 
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Figure 3: (Color online) Non-interacting model with U a = Ud = 0, 
asymmetric case, 7^ 0. a) Conductance curves for a range of tem- 
peratures. Full curves: NRG results. Dashed lines: Fano resonance 
line fits, b) Fano parameters as a function of the temperature. Pa- 
rameters are extracted by curve fitting in the energy window which 
corresponds to the horizontal axis in the upper subfigure. 



(by one electron for each spin orientation) as e a crosses the 
Fermi level, thus the phase change in each spin channel is tt: 
this corresponds to the usual Fano-resonance scenario. In the 
second limit, the total occupancy changes only by one due 
to the electron-electron repulsion preventing the second elec- 
tron from entering the quantum dot. This corresponds to a 
7r/2 phase shift in each spin channel and, roughly speaking, 
only one half of the Fano-resonance-like feature is expected 
to be seen. In this section we discuss the first case, i.e., the 
U a *C r a limit, which we study by setting U a to zero. 

When the embedded dot is tuned to its particle-hole sym- 
metric point (Sd = 0), the impurity spectral function A° d is 
pinned at the Fermi level to the value of 1/ 7rr d irrespective 
of the value of the electron-electron repulsion Ud- This is true 
even for moderate departures from the particle-hole symmet- 
ric point, \8d\ < Ud/2, since the Kondo resonance remains 
near the Fermi level, unlike the Lorentzian peak in the non- 
interacting case. For small t, the characteristic energy scale of 
the side-coupled quantum dot is therefore still T a = t 2 /Td 
as in the symmetric non-interacting case. For larger t the 
broadening will, however, eventually become comparable to 
the width of the Kondo resonance (of order Tj ( ). The behav- 
ior therefore depends on the ratio between the hybridization 
of the side-coupled dot r a and the Kondo temperature of the 
decoupled embedded dot, T%. If T a < T^, the physics of the 
embedded dot will not be affected much by the side-coupled 
dot (and vice versa), thus the conductance as a function of e a 
will differ only slightly from the non-interacting case. This is 
illustrated in Fig. [4] upper panel, where we compare conduc- 
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tance curves for U d = and U d /D = 1 at f/D = 0.0001, 
where r a /T^- = 1.7 x 10~ 3 since the Kondo temperature 
(Wilson's definition) is T% = 9.7 x 10~ 5 L>. In this regime 
the Kondo effect on the embedded dot in no way affects the 
Fano interference process in the relevant range of energies and 
temperatures (i.e. several times r a ). Only for very large e a 
and very high temperature (of the order of T^) will the differ- 
ences become apparent, however this is outside the parameter 
regime of interest here. 

When T a is equal to a considerable fraction of T^-, we start 
to see small quantitative departure from the Fano line shape 
already for e a and T of the order of T a , see Fig. @] middle 
panel. This is clearly a consequence of the competition be- 
tween correlation and interference effects, combined with fur- 
ther thermal effects. As expected, the discrepancy grows with 
increasing e a and T. 

Finally, for T a 3> T%, the embedded dot is strongly per- 
turbed by the coupling to the side-dot and the differences 
become drastic: at finite T, the line-shape differs qualita- 
tively from the Fano form and, in particular, we observe the 
emergence of non-uniform independence with broad humps 
around E ~ T„ . 
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Figure 4: (Color online) Resonance curves for non-interacting (full 
lines) and interacting models (dashed lines) for a range of scaled tem- 
peratures T/r a and for different Ta/Tp* ratios. 

It is instructive to perform curve fitting on the results where 
the Fano resonance is already strongly perturbed, as shown 
in Fig. The fits (shown using dashed lines), performed 
in the energy interval [— 4r a : +4T a ], are clearly inade- 
quate. Extracted Fano parameters indicate that the resonance 
width is significantly reduced below r a and its width even 
decreases with increasing temperature, which reflects the sit- 



uation where the spectral function of the embedded dot has a 
significant variation on the energy scale of r a , which affects 
the resonance line-shape for large E. It should also be noted 
that parameters a and b have "unphysical" values a > 1 and 
b < for small T. If the curve fitting is performed in a nar- 
rower energy interval [— T a : +T a ], the asymptotic small-i? 
form of the resonance can be well captured, however the ex- 
tracted Fano parameters are factitious and therefore of little 
use. This demonstrates that in the presence of strong com- 
petition between the Fano interference and Kondo effect, the 
curve fitting using the Fano line-shape is not recommended 
since the results depend very strongly on the energy window 
where fitting is performed and it is thus not advisable to make 
any inference based on them. 
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Figure 5: (Color online) Interacting model with Ud ^ and U a = 
0. Symmetric case, ed = 0. a) Conductance curves for a range of 
temperatures, b) Fano parameters as a function of the temperature. 
Curve fitting is performed in the energy window which corresponds 
to the horizontal axis in the upper subfigure. 

Since there are two different energy scales in the problem 
{Tj ( and T a ), the temperature-dependence of the conductance 
is expected to be non-monotonic, as shown in Fig.|6]where we 
plot the conductance at the bottom of the Fano anti-resonance 
(i.e. for e a — 0). Such dependence is a consequence of the 
competition between the Kondo and Fano effects. The Kondo 
effect tends to increase the conductance through the formation 
of many-particle resonance at the Fermi level which opens a 
new conduction channel through the system. On the other 
hand, the Fano effect suppresses the conductance through 
quantum interference. For Tjj ;§> T a , the conductance first 
increases at the higher temperature scale of Tjj (the tempera- 
ture dependence being given by the universal Kondo conduc- 
tance curve 71 ), then it decreases at the lower temperature scale 
set by T a : see results for t/D = 0.0001 and t/D = 0.001 in 
Fig. [6] Note that the unitary limit of full conductance quantum 



6 



is not achieved at intermediate temperatures even in the case 
of well separated energy scales (three orders of magnitude for 
t/D = 0.0001). For larger t, the conductance peaks near 
T ~ r a , where both competing effects are equally suppressed 
by thermal fluctuations. 



1 




T/T 

a 

Figure 6: (Color online) Temperature-dependence of the conduc- 
tance for e a = 0, i.e. at the bottom of the conductance anti- 
resonance, for a range of the inter-dot coupling strengths. 

It is also worthwhile to compare the temperature- 
dependence of the conductance G(T) with the energy- 
dependence of the zero-temperature spectral function 
Ad(w,T = 0), see Fig. [7] A common approximation in 
discussing the transport properties of quantum dot systems is 
to assume that the temperature variation of the conductance 
simply follows the energy variation of the zero-temperature 
spectral function, i.e., G{T) ra nY d A d (T, 0). While this 
is true on the qualitative level, the numerical results shown 
significant qualitative differences: the conductance curves 
are generally significantly broader and, accordingly, their 
maxima are lower. In particular, while the spectral function 
often attains values very close to 1/nTd in some frequency 
range, the conductance is not unitary in the corresponding 
temperature range. 

The parameter regime studied in this section corresponds 
to the model used in Ref. 34 to obtain the conductance plots 
featuring an anti-resonance, which were proposed to explain 
the experimentally observed Fano-resonance-like line shapes 
with small q. It can be noticed, however, that this parameter 
regime is not appropriate, since it corresponds to a change of 
occupancy of the side dot by 2, while from the experimentally 
observed honeycomb stability diagram one can deduce that 
the electron number changes only by one. We return to this 
issue in Sec. [VI] 



V. FULLY INTERACTING CASE 

We first study the effect of finite electron-electron repul- 
sion on the side-coupled dot in the case where the embed- 
ded dot itself is non-interacting, U d = 0. In this way we 
exclude all processes related to the inter-dot exchange inter- 
action and we may focus on the path-interference physics. 




Figure 7: (Color online) Comparison of the temperature-dependence 
of the conductance G(T) with the energy-dependence of the zero- 
temperature spectral function on the embedded dot Ad(w,T — 0). 
For completeness, the spectral function on the side-coupled dot 
A a (u),T = 0) is also shown. Spectral functions are rescaled in 
units of the respective effective hybridisation strength, i.e., by l/7rri, 
where i £ {a, d}. Conductance curves are rescaled in units of the 
conductance quantum Go = 2e 2 /h. The horizontal axis is rescaled 
in units of F a — t /Yd (note that r a is different in each subfigure). 



To simplify the discussion, we focus solely on the symmet- 
ric case with = 0. As previously discussed, as long as 
U a <C r a , the side-coupled dot is essentially non-interacting 
and the results are equivalent to those shown in Sec. lIII Al For 
larger U a > irT a , however, the side-coupled dot will undergo 
the Kondo effect and its occupancy will be pinned to 1 for e a 
below several times r a . Consequently, the "resonance" line- 
shape in this regime will take the form of an inverted "Kondo 
plateau onset". The evolution from the Fano resonance to the 
Kondo plateau onset is shown in the upper left panel in Fig. [8] 
while the associated level occupancy curves are shown in the 
bottom left panel. (See also Fig. 3 in Ref. 53 where a similar 
setup is considered.) 

We now finally discuss the fully interacting case, which is 
the situation which is relevant for the actual experimental con- 
figuration. The results are shown in the right panels in Fig. [8] 
We find a very similar evolution in the conductance and occu- 
pancy curves with increasing U a as in the Ud = case. The 
transition from the Fano line shape to the Kondo plateau-edge 
shape occurs in this case for lower values U a , however, and al- 
ready by U a /D = 0.01 we obtain conductance curves which 
do not change much as U a is further increased. This is very 
important since it implies that the regime of significantly large 
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Figure 8: (Color online) Conductance curves and level-occupancies 
at T = for a range of electron-electron interaction strengths U a ■ 
Left panels: non-interacting embedded impurity. Right panels: in- 
teracting embedded impurity. The arrows show the direction of in- 
creasing U a - 

U a occurs for rather small values of U a , therefore it is not ap- 
propriate to model the side-coupled dot as a non-interacting 
impurity, since it is unlikely that the electron-electron repul- 
sion in two similar quantum dots would differ by orders of 
the magnitude. We now face a dilemma: assuming a non- 
interacting side-coupled dot, we obtain conductance curves 
which appear to agree with the experimental results (albeit we 
know that they correspond to a change of occupancy by 2, in 
contradiction with the experimental situation), while assum- 
ing an interacting side-coupled dot, we obtain very different 
results which are in strong disagreement with those actually 
observed. 
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Figure 9: (Color online) Comparison of the temperature-dependence 
of the conductance G(T) with the energy-dependence of the zero- 
temperature spectral function on the embedded dot Ad(u}, T = 0). 

In Fig. [9] we compare the temperature dependence G(T) 
with the energy-dependence of the T = spectral function 
A d (uj,T = 0). As before (cf. Fig. Q and the related dis- 
cussion), we observe that in spite of a certain general simi- 
larity between the curves, the differences are notable. In this 



case, the differences are even qualitative: while G(T) exhibits 
two ranges of logarithmic behavior (the first Kondo screening 
stage from 10~ 4 < T < 10 ~ 2 , and the second Kondo screen- 
ing stage from 10 -6 < T < 10~ 4 ), the zero-temperature 
spectral function is more complex and even attains a minimum 
at uj « 2 x 10~ 6 before increasing back to its asymptotic value 
at the Fermi level. This complex behavior is not reflected in 
the temperature dependence of the transport properties. Ad- 
ditional calculations (results not shown) indicate that G(T) 
and A^(uj,T = 0) are very similar only in vicinity of the 
particle-hole symmetric point, and even there only when the 
two Kondo temperature scales are well separated. This should 
serve as a caveat: since G(T) is a convolution of the finite- 
temperature spectral function, see Eq. (0, it is not expected 
that, in general, all the details of the energy-dependence of 
the zero-temperature spectral function will manifest as similar 
variations in the temperature-dependence of the conductance 
at the related temperatures. 

VI. RELATION TO THE EXPERIMENTAL RESULTS 

Since the experiment is performed at finite temperatures, 
the results for the zero-temperature conductance can be 
very misleading in problems with very low energy scales 50 . 
This commonly occurs in the side-coupling geometry and 
in clusters of quantum dot with weak inter-dot tunnel 
couplin g 27 ' 28 ' 48 ' 50 . In Fig. [TUh ) we thus plot the conductance 
for a range of finite temperatures, which reveals completely 
different transport properties when compared with the zero- 
temperature limit. 

When both quantum dots are interacting and the parameters 
are set in such a way that local moments form on both dots, 
there is a possibility that the Kondo screening will proceed 
in two step a 26,28 ' 29 ' 48 ' 50 ' 52 ' 73 : at higher Kondo temperature T^, 
the moment on the first dot will be screened, while another 
stage of Kondo screening of the moment on the side-coupled 
second dot occurs at another much lower energy scale Tj p 8 ' 73 : 

T 2 K = ciT° K exp (-c 2 T° / J cfi ) , (15) 

where c\ and c 2 are some constants of the order 1 . Such two- 
stage Kondo screening occurs if the effective Kondo exchange 
coupling Jeff of the spin on the second dot with the quasipar- 
ticles resulting from the first stage is lower than In this 
case, the second Kondo scale may be pushed to exponentially 
low temperature Tf ( . The exchange coupling 

Jeff = -JT- + -JT- (16) 

e Q - £d + Ua e d - e a + Ud 

depends on the value of e a and decreases as e a is reduced, 
attaining its minimum for e a = 5d — U a /2, In turn, this im- 
plies that becomes very low for small e a , thus at the ex- 
perimentally relevant temperatures the second stage of Kondo 
screening does not occur, the conductance suppression does 
not arise and the conductance curve takes a form similar to 
that of the Fano resonance. In experiments, the temperature 
is roughly one order of the magnitude lower than the Kondo 
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Figure 10: (Color online) a) Conductance curves for the fully inter- 
acting case with Ud/D — U a /D = 1, 5d = 0, for a range of differ- 
ent temperatures, b) Temperature dependence of the conductance for 
a range of energies of the side-coupled dot. 



temperature; for chosen t/D = 0.001, we indeed find that for 
T/D = 10" 5 « O.lTg-, the conductance form takes the form 
of a nearly symmetric anti-resonance, in agreement with the 
experiment. The minimum of the anti-resonance is located at 
lower value of E = e a /T a as compared to the non-interacting 
case (where it is found exactly at E = 0), since the physi- 
cal mechanism for its occurrence is completely different: in 
the non-interacting case we observe a complete suppression 
of the conductance due to quantum interference between dif- 
ferent conductance pathways, while in the interacting case the 
incomplete suppression results from the separation of the en- 
ergy scales Tj{ and with the experimental temperature set 
on an intermediate scale so that < T < T^. The shift 
from E = is not observable experimentally, since the exact 
value of the parameter e a is not known: it has to be inferred 
indirectly from known gate voltages on the electrodes and the 
capacitances. 

If the inter-dot tunneling t is increased so that the exchange 
coupling becomes comparable or exceeds the Kondo temper- 
ature, the conductance curves no longer exhibit a Fano-like 
antiresonance, not even at elevated temperatures, see Fig.fTTI 
Instead, as E is reduced so that the side-coupled dot becomes 
occupied by one electron, the two electrons (one from each 
quantum dot) bind to a local singlet on the temperature scale 
of J ~ At 2 1 U and there are no other lower temperature scales. 

The interpretation of the anti-resonance in terms of a ther- 
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Figure 1 1 : (Color online) Conductance curves for the fully interact- 
ing case with Ud/D — U a /D = 1, 5d = 0, for increasing inter-dot 
tunneling t. The arrow indicates the direction of the increasing tem- 
perature T. 



mally suppressed second-stage Kondo screening is supported 
by the experimental results for the temperature-dependence of 
the resonant conductance (Fig. 2c in Ref. 34) which indicate a 
logarithmic behavior. This temperature variation corresponds 
to the high-temperature range in Fig. [T0b). The experimen- 
tal results do not indicate any saturation of G for the lowest 
experimentally accessible temperatures, however this might 
simply indicate that the temperature is not sufficiently below 
Tk, or that t is indeed very small and thus T\ is exponentially 
low. Further confirmation for the interpretation in terms of 
the two-stage Kondo effect could be obtained by also consid- 
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ering the temperature-dependence of other parameters which 
determine the anti-resonance line-shape; the variations of am- 
plitude, width, and asymmetry parameter should all exhibit 
logarithmic temperature dependence. 

In addition to the thermal effects, the second stage of Kondo 
screening may also be suppressed by a minute magnetic field 
such that I|<B<3|, As it is likely that there are al- 
ways some small stray magnetic fields in the experimental de- 
vice, this provides another mechanism which tends to produce 
Fano-resonance-like spectral features in the present configura- 
tion. 



VII. CONCLUSION 

The transport properties of a system of two quantum dots 
in the side-coupled geometry have been analyzed, focusing on 
the parameter range where the occupancy of the side-coupled 
dot changes as its energy crosses the Fermi level. The re- 



sulting resonance line-shapes in the conductance curves have 
been comprehensively studied in various parameter regimes 
for both non-interacting and interacting quantum dots. We 
find an alternative interpretation of the experimental results 
from Ref. [34j which fulfils the two required conditions: 1) 
the occupancy in the side-coupled dot changes by one, 2) the 
resonance line-shape features a weakly asymmetric Fano-like 
anti-resonance. The proposed model takes into account the 
electron-electron interaction on the side-coupled dot, which 
is certainly present in the experimental device and which has 
been shown in this work to be relevant for the transport prop- 
erties of the system. 
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